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A calculation of the QCD phase diagram at finite temperature, 
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We study the phases of a two-flavor Nambu-Jona-Lasinio model at finite tempera- 
^ \ ture T, baryon and isospin chemical potentials: fis = (/^M+Md)/2 , fii = (fi u — fj,d)/2. 

This study completes a previous analysis where only small isospin chemical potentials 



Hi were considered. 
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The possible formation of a pion condensate due to a finite isospin chemical "potential 
P-i! jx^ has been, in recent years, the subject of several papers |l|, 



12 



13. 



14J. Consequently, the reconstruction of the QCD phase diagram at finite tempera- 
ture and quark densities, such as those attainable in earth experiments and in the interior 
of stars, is even more challenging due to the addition of /i/, besides temperature T and 
baryon chemical potential hb- Various regions of the phase diagram correspond to different 
experimental settings. Actually, the behavior of QCD at high temperature and low baryon 
densities is central to the relativistic heavy-ion collisions: experiments at CERN and RHIC 
are expected to produce hadronic matter in this regime. On the other hand, the description 
of neutron star interiors requires the knowledge of cold nuclear matter at large baryon den- 
sities. However, nature does also provide us with systems at finite isospin chemical potential 
Hi in the form of asymmetric-isospin matter inside neutron stars: nuclear matter has a finite 
(negative) isospin 7 3 density due to Coulomb interactions, apart from finite baryon-number 
density. Moreover, in any realistic experimental setting in relativistic heavy-ion collisions 
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there is a non-zero, even if small, 

Our present description of the QCD phase diagram in the plane (//b,T) anticipates the 
existence of a tricritical point separating first order transitions in the regions of low temper- 
atures from cross-over transitions in the low baryon chemical potential and high temperature 
regime [3 [3 12, 13 13 13 13] • In recent years various other non-trivial phases of QCD 
at low temperatures and high baryon chemical potentials have been discovered, such as the 
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jocking (CFL) phase and the two-color superconducting (2SC) phase (for a 
)• 

Coming back to the effects of the isospin chemical potential fii, pion condensation has so 
: ar been primarily investigated by means of low energy models based on chiral lagrangians 
sl, Q, 0, I3|- Although these models are well suited to study the phases of QCD as they 
have the right symmetry properties, they do not include the combined effects of the isospin 
chemical potential fij with a finite baryon chemical potential /i# in order to study the pattern 
of chiral symmetry breaking and restoration as well. To consider both fiB and /i/, we need 
a model with quarks as microscopic degrees of freedom. The effect of small fii (up to half of 
the pion mass) has been investigated in ref. 3| in the context of the Nambu-Jona-Lasinio 
model and in ladder-QCD 14]. The result is the splitting of the critical curves for chiral 
symmetry restoration for the two light flavors, whereas a full study for arbitrary fij has only 



been done in the context of a random matrix model 11 U . 

Studies on the lattice have been performed at finit e uj and an = in refs. |3 25, 13 
27|, 3] and with a finite /ib and /if = in refs. jssL [20!. I31L [3^] . In a recent work 3| the 
effect of both \xb and a small \ii has also been considered. The case of high fiB and small 
Hi has been considered in (34 1 . 

In this work we extend the analysis of HQ where it was found that the first order 
transition line ending at the tricritical point of the case fij = actually splits into two first 
order transition lines and correspondingly two crossover regions are present at low values 
of baryon chemical potential. In particular we will be working in the context of a NJL 
model with a form factor included such as to imply a decreasing of the fermion self-energy 
compatible with the operator product expansion. 

It should also be noticed that in [3] the NJL model has been augmented by the four- 
fermi instanton interaction relevant in the case of two flavors. These authors have found that 
the coupling induced by the instanton interaction between the two flavors might completely 
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wash the splitting of the first order transition line. This happens for values of the ratio of 
the instanton coupling to the NJL coupling of order 0.1-0.15. 

In Section II we summarize the relevant features of the NJL model we have considered, 
with isospin charge included. The one-loop effective potential and the values of the fit param- 
eters are included. In Section III we discuss the various equilibrium phases of the model, 
together with the corresponding symmetries, by studying the behavior of the scalar and 
pion condensates with respect to different thermodynamical parameters among T, /xg, fij (or 
fi u , fid). Results are shown for growing temperatures, starting from zero up to temperatures 
above that of the critical ending point. Finally, Section IV is devoted to conclusions. 



II. THE MODEL 



Our purpose is to explore the structure of the phase diagram for chiral symmetry and 
pion condensation in QCD at finite temperature and quark densities, by using a microscopic 
model with quark degrees of freedom. This task has been accomplished, up to now, in the 
context of a random matrix model simulating QCD with two flavors Q and, in the case of 
small differences between the u and d quark chemical potentials, also in the Nambu-Jona- 
Lasinio model (NJL) and in ladder-QCD 1^ . 

One reason for using a model with quarks as microscopic degrees of freedom is that it gives 
us the possibility of studying chiral symmetry breaking and pion condensation at both finite 
isospin and baryon chemical potentials, which is not possible within effective chiral models. 



In ref. |12[, the authors made use of the NJL model with a suitable form factor included 
in the quark self-energy to mimic asymptotic freedom 0, 0]. This version of the NJL 
model turns out to be very close to ladder-QCD as developed in refs. where the 

momentum dependence of the quark self-energy is consistently dictated by the study of the 
Schwinger-Dyson equation within a variational approach (see the previous references for 
details). However, although ladder-QCD is a covariant and self-consistent approach, the 
dependence on the four-momentum of the quark self-energy makes the numerical computa- 
tion of the one-loop effective potential with finite quark densities much more onerous with 
respect to the NJL case, where the quark self-energies depend only on the three-momentum. 
For this reason, in the present work we study the NJL model. It is reasonable to expect 
that when employing ladder-QCD , the resulting physical picture does not considerably differ 



from that of the NJL model. This has been the case in previous applications too 16|,ll7|, 3|. 



As already said in the Introduction we are not going to consider the effects of di-fermion 
condensation. Therefore our results can be considered valid only outside the region of the 
color superconductive phase, which is roughly in the region defined by fis ^ 400 — 500 MeV 
and T < 50 MeV. At the same time we will not consider regions at values of fi u or Hd higher 
ihan 400 — 500 MeV where other difermion condensates might arise (see for instance ref. 



m 

Let us now consider the Lagrangian of the NJL model with two flavors u, d with the same 
mass m but different chemical potentials /i u and fid 

£> — £o + Cm + + Li n t 
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tyidty - m W + ^ f A * + — ^2 [{^ r a^) 2 + (^75^) 2 (1) 





G 

a=0 



where \1/ = I , A = I is the matrix of chemical potentials and 

\d J \0 lid) 

r a , a = 0, 1, 2, 3, is the set of the three Pauli matrices plus the identity. 

We note that we can express £ M either by using the variables fi u , fid or the two combinations 

Hb = — ~ ^ d an d fii = ^" - ^ d , which couple to the baryon charge density and to the 

third component of isospin respectively 

£ M = fi B tfttf + m ^t T3 ^ ( 2 ) 

To study whether a pion condensate shows up, we need to calculate the effective potential. 
This is obtained by using the standard technique to introduce bosonic (collective) variables 
through the Hubbard-Stratonovich transformation and by integrating out the fermion fields 
in the generating functional. However, the effective potential that we have considered is not 
directly obtained from the Lagrangian in Eq. (0). To mimic asymptotic freedom we want 



;o include a form factor as in ref. |35l and we thus follow the same procedure as in refs. 

in 

12L 12 II ] . The result is a one- loop effective potential which generalizes that of the theory 
described by the Lagrangian in Eq. (|I|L and which reduces to it in the limit of a constant 
form factor F(p) = 1. 

V=^(xl + X 2 d + 2p 2 ) + V lQg (3) 
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h u -F 2 (p) A p 7 5 
F 2 (p) A p 7 5 h d 

(4) 

h f = (iu n + y u / )7 - p-7 - (m + F 2 (p) A Xf) 

where c<j n are the Matsubara frequencies and the dimensionless fields Xf an d P are connected 
to the condensates by the following relations 

Xf = -2G <*'*/> 



A 

(5) 

(tt7 5 (i - ^7 5 m) 
A 



and are variationally determined at the absolute minimum of the effective potential. In the 

A 2 fl 

previous equations, A is a mass scale appearing in the form factor F(p 2 ) = ^ 2 ^ — - |35j . 
It is worth noting that the one-loop effective potential in Eq. (jlj) has the same expression 
of the one derived in ref. j3| within ladder-QCD. Therein, multiplying the scalar and 
pseudoscalar fields, there was a test function guessed from the study of the one-loop 
Schwinger-Dyson equation for the quark self-energy, in place of F 2 (p) in Eq. (j3J). The only 
difference is that F 2 depends on the three-momentum whereas the quoted test function 
depends on the four-momentum and that the two asymptotic behaviors are different 
(~ 1/p 2 in the test function of ref. ^| and ~ 1/p 4 in F 2 of Eq. (0])). Otherwise the 
two effective potentials would be identical. This observation also explains the reason why 
we have adopted the NJL model instead of ladder-QCD to generalize the analysis of ref. 
Q| at high isospin chemical potentials as the numerical analysis is much simpler in this case. 

To fix the free parameters of the model, which are A, the average current quarks mass 
m = (m u + md)/2 and the coupling G = g/A 2 , we work at zero temperature and quark 
densities. We first choose the mass scale A within the range A ~ 500 — 600 MeV. Then 
we determine the strength of the coupling g and the mass parameter m by requiring a light 
quark condensate of the order (^ff^ff) ~ —(200 MeV) 3 and a pion mass ~ 140 MeV 
(the latter evaluated through the curvature of the effective potential in the direction of the 

n 

pion field and having fixed f n at its experimental value [36]). 
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The output parameters are the following 

A = 580 MeV; g = 7; m = 4.5 MeV (6) 

With these values we obtain a condensat e (J &ji&f) = —(172 MeV) 3 and a constituent 



quark mass Mf = 428 MeV (defined as in |l9j,|2l(). The critical isospin chemical potential 
at zero temperature turns out to be pf = 89 MeV. The discrepancy of about 25% (we 
recall that the expected value of pf would be m n /2 ~ 70 MeV) is due to the approximate 
fit procedure. Actually, in our previous work based on ladder-QCD 14], where f n was 
consistently calculated within the model, we got pf = m 7T /2. 



III. PHASE DIAGRAM FOR CHIRAL SYMMETRY BREAKING AND PION 

CONDENSATION 

In order to discuss the structure of the phase diagram, it is worth summarizing the 
symmetries of the Lagrangian density in Eq. ((TJ). Both C and Ci nt are SUl(2) <S> SUr(2) ® 
Uy (1) Cg) U^(l) invariant. The symmetry is reduced by the mass term C m to SUv{2)®Uy (1) 
and further reduced from the term £ M which selects a direction in the isospin space, as is 
evident from Eq. (j2J), unless p u = pd and thus fii = 0. The remaining symmetry can be 
expressed either as Uy(l) <g> U V {1) or Uy(l) <g> U v (l), depending on the basis of the fields 
that we are choosing. 

The baryon number symmetry Uy(l) is dynamically respected, whereas a non vanishing 
v.e.v. of the p field defined in Eq. (jHJ) may appear, which dynamically breaks U v (l). This 
implies the appearance of a Goldstone mode, which is either the charged tt + or at the 
threshold, depending on the sign of pi, whereas the other two pions are massive modes. 
As far as the scalar condensates Xu,Xd (see Eq. ©) are concerned, they do not break any 
symmetry. However, since the mass term is small, their value is almost entirely due to the 
approximate spontaneous breaking of chiral symmetry. Consequently we distinguish regions 
where the dynamical effect is relevant, from regions where the scalar condensates are of order 
~ m/A, namely where only the effect of the explicit breaking of chiral symmetry survives. 
The determination of the various phases has been performed numerically by minimizing the 
one-loop effective potential. We start by showing the results in the (p u , Pd) plane, for fixed 
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values of the temperature. Different regions are labelled, as in ref. ll|, by the symbol of 
the field which acquires a non vanishing v.e.v. due to dynamical effects, whereas the other 
fields are vanishing (p), or of the order ~ m/A (x u and/or Xd)- 

Solid lines refer to discontinuous transitions and dashed lines to continuous ones. However, 
we recall that strictly speaking only the lines surrounding regions with a non vanishing field 
p refer to genuine phase transitions, associated with the breaking and restoration of the 
Uy(l) symmetry. 



A. Zero temperature 

In Fig. ^ we show the phase diagram in the (p u , Pd) plane at zero temperature. Let 
us start from the vacuum at T = p u = pd = 0, at the center of the picture. Here, in the 
chiral limit, the pions are the Goldstone bosons associated with the spontaneous breaking 
of SU{2)l ® SU(2)r. We have chosen these variables in order to compare the structure of 
the phase diagram with that obtained in ref. However, if we want to recover known 

results given in terms of the baryon chemical potential, we have to move from the center 
along the diagonal at p u = pd and thus at p I = 0, and increase the absolute value of ps- 
At — ^ ^ = \fig\ = 293 MeV we meet the approximate restoration of chiral symmetry 
due to the sudden jump of the condensates of the two (degenerate) quarks to values of order 
~ m/A, which is a discontinuous transition. The same thing happens by moving along lines 
parallel to the main diagonal in the region labelled by Xu, Xd, enclosed between the two 
dashed lines at — - ^ = \pj\ = 89 MeV (see also Fig. El where it is shown that the 
critical value of pj at T = is independent on Pb) and by varying \pb\- The regions in 
the top-right and bottom-left corners of Fig. ^thus have the U v (l) <g> U v (l) symmetry of £ 
with p = and Xu,Xd of order ~ m/A. 

By moving from the center along the diagonal at p u = —pd (and thus ps = 0) or parallel to 
it, when crossing one of the two dashed lines at \pi\ = 89 MeV (we have already discussed 
the origin of the discrepancy between this value and half of the pion mass in the model), the 
absolute minimum of the effective potential starts to rotate along the p direction. We thus 
have a continuous breaking of U v (l) and a second order phase transition with one Goldstone 
mode which is, right along the dashed line, either the ir + (in the upper part of the diagram) 
or the 7r~ (in the lower part). In the chiral limit these two dashed lines merge together 
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in coincidence with the diagonal at pi = as the pion becomes massless in this limit and 
the rotation is sudden, giving first order phase transitions for pion condensation. In this 
case there are two Goldstone bosons associated with the spontaneous breaking of two U(l) 

. 

Coming back to the massive case, and still with reference to Fig. ^ we conclude that by 
considering \ps\ not too large and by growing \pi\, we find a second order phase transition 
with the rotation of the scalar condensates into the pseudoscalar, namely we are faced with 
pion condensation in a 

not find the vanishing of p for values of \pi\ high with respect to the pion mass, but still 
sufficiently low to avoid considering superconductive phases (actually, for very low \ps\ this 
transition would occur in the present model for \pi\ ~ 1 GeV). 

To explore the possibility of multiple phase transitions and thus of a richer phenomenology, 
we need to grow \pb\ as for instance we do in the case described in Fig. |3]where we follow the 
path of the solid line a in Fig. [T]at ps = 170 MeV for growing pi > 0. The fields Xu and Xd 
are almost degenerate, both in the region of the approximate dynamical breaking of chiral 
symmetry (below pj = 89 MeV) and in the region of spontaneous breaking of U v (l), where 
they rotate into the p field. Then, when the line a in Fig. ^crosses the solid line surrounding 
the region labelled by Xd, we see that p suddenly jumps to zero with the restoration of the 
Uy(l) group and that Xu and Xd split. Actually the latter suddenly acquires a value due 
to the dynamical breaking of chiral symmetry whereas Xu undergoes a further decrease and 
remains of order ~ m/A. 

In Fig. |Uwe plot the behavior of the scalar condensates Xu and Xd vs. pi at ps = 210 MeV, 
namely by following the path described by the solid line b in Fig. ^ We see that we never 
cross the region with p ^ and that we simply pass from a region where the dynamical 
effect of the breaking of chiral symmetry is entirely due to a large value of Xu and Xd of 
order ~ m/A (at large negative pi and small p u ), to a region where this effect manifests 
itself with a large value of Xd and Xu of order ~ m/A (at large positive pi and small pd). 
The region in between has almost degenerate and both large Xu and Xd- Pion condensation 
does not occur for this value of ps (see also Fig. EJ). 

Finally, in Fig. |3J we plot the behavior of the condensates at fixed p u = 200 MeV vs. pd 
(see again Fig. The rotation of the pion condensate into the scalar ones occurs when the 
vertical line at p u = 200 MeV meets the dashed line at pj = 89 MeV, which happens for 
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Pd of few MeV . Then, when pd has sufficiently increased, Xd falls to a small value of order 
~ m/A with a discontinuous transition, whereas Xu remains constant at its large value. 

B. Finite temperature 

The evolution of the phase diagram for growing temperatures is easily understood as far 
as the regions with p = are concerned. Actually, in this case the effective potential at 
the minimum is the sum of two independent terms, one for each flavor, and the results are 
straightforwardly given through the analysis of chiral symmetry breaking and restoration 
or a single flavor at finite temperature and chemical potential (see for instance refs. 



Ill (3 )- In Fi g. El we show the phase diagram at zero, or small isospin chemical potential 



lg. |Qjwe f 

HQ; 



(see also refs. |12L 1141]) . From this picture we see that moving along any of the critical 
lines of chiral symmetry restoration at fixed pi, the critical value of the baryon chemical 
potential ps decreases for growing temperatures. Furthermore, for temperatures below 
that of the critical ending point E, T < T(E) = 85 MeV, the transitions are always 
discontinuous whereas they become cross-over transitions for T > T(E). Consequently the 
regions labelled by Xu and/or x<i in Fig. ^ shrink when growing T and their rectilinear 
sides become lines of cross-over transitions for T > T(E) (see Fig. Hand Figs. |H1 where 
we plot the phase diagrams in the (p u , pa) plane at T = 60 MeV and T = 140 MeV, which 
is respectively below and above T(E)). The new feature concerns the regions with p ^ 0, 
which also reduce their size for growing T, whereas the order of the transitions starts to 
change from first to second, beginning from the critical points at highest \pi\, until they 
reach the points of the boundaries which coincide with those of the regions labelled by 
Xu or Xd (see Fig. |HJ). Also the length of the curves of second order phase transitions to 
pion condensation at fixed values of \pb\ sensibly reduces from low temperatures to high 
temperatures (see again Fig. and Figs. 0IHJ). A similar behavior, for high T, is found in 
ref. nj. For T > T(Pj-) = 174 MeV, which is the cross-over temperature at zero chemical 
potential (see Fig. EJ, all these regions disappear from the phase diagram, which is thus 
characterized by p = and Xm Xd ~ m/A. 

The behavior of the scalar and pseudoscalar condensates at T = 60 MeV are much similar 
to those at T = 0. As an example we plot, in Fig. El the condensates at p u = 200 MeV 
and T = 60 MeV vs. pd (compare with Fig. EJ. The situation is different if we consider 
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temperatures above T(E) = 85 MeV. For instance, at T = 140 MeV, we see from Fig. 
IH1 that the structure of the phase diagram is only slightly modified with respect to the 
case of two independent flavors which undergo cross-over phase transitions at sufficiently 
high values of their own chemical potentials (actually the region of pion condensation has 
sensibly reduced with respect to Figs. ^ HJ). The phase transition associated with the 
spontaneous breaking of Uy(l) can be both second or first order, depending on the path 
followed. In Fig. |H1 the solid line a refers to a path at fig = vs. \ii > where the 
transition to pion condensation is continuous. The behavior of the condensates relative to 
this path is plotted in Fig.. EH In Fig. [TT]we plot the scalar and pseudoscalar condensates 
at \i u = 200 MeV and T = 140 MeV vs. fi d . 

In Fig. El we plot the value of the critical isospin chemical potential fif beyond which a 
pion condensate forms vs. temperature T at zero baryon chemical potential fiB- The growth 
of \xf is easily understood since the pion mass (which should be twice [if ) is expected to 
grow near the critical temperature for chiral symmetry restoration On the other hand, 
no phase transition to pion condensation is expected above the cross-over temperature for 
chiral symmetry restoration. Thus the line of critical values ends at T = 174 MeV. 
The pion condensate is also expected to decrease for growing temperatures: in Fig. [T31 we 
show p vs. T at /is = and fii = 200 MeV. Similar behaviors are obtained for fixed values 
of fiB and fij. 



IV. CONCLUSIONS 

In this paper we have continued the study of pion condensation at finite baryon and 
isospin density in the NJL model that we started in ref. Q| in the case of ladder-QCD 
for small isospin chemical potentials. The extension to higher isospin chemical potentials 
confirms the structure predicted in ref. n|], where two-flavor QCD was simulated in the 
context of a random matrix model. Some difference between the two analysis is present, 
at low temperatures, in the region of high isospin chemical potentials, at the boundary 
of the region where color superconductivity should take place. Actually in this region we 
find that pion condensation is still active, whereas in ref. [ljj the authors find that the 
pion condensate vanishes. We have also shown the expected behavior of scalar and pion 
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condensates by following different paths, for growing temperatures, both in the plane of 
quark chemical potentials //^) and in that of isospin and baryon chemical potentials 
(///, ijlb)- The analysis that we have performed should also be confirmed, with only small 
quantitative differences, within ladder-QCD . 
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FIG. 1: Phase diagram for chiral symmetry restoration in the plane [puil^d] of quark chemical 
potentials, at T = 0. Different regions are specified by the non vanishing of a given condensate, 
whereas the others are vanishing (p) or order ~ m/A (\ u an d Xd)- Dashed lines are for the 
continuous vanishing of p or for cross-over phase transitions for Xu or Xd> whereas solid lines are 
for discontinuous behaviors. The solid lines a and b refer to specific paths at fixed values of ps, 
with pb = 170 MeV (line a) relative to Fig. \^and ps = 210 MeV (line b) relative to Fig. ^[ 
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FIG. 2: Critical value of the isospin chemical potential, beyond which a pseudoscalar condensate 
forms, vs. baryon chemical potential, at zero temperature, [ib = 204 MeV is the highest allowed 
value for pion condensation to occur. The path followed in the phase diagram of Fig. Q] is along 
the upper-half of the dashed line in the lower half-plane. 
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FIG. 3: Scalar and pseudoscalar condensates vs. m, for ^b = 170 MeV and T = 0. The path 
followed in the phase diagram of Fig. Q is that of the solid line a. 
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FIG. 4: Scalar condensates vs. [ii for \jlb = 210 MeV and T = 0. The pseudoscalar condensate is 
zero. The path followed in the phase diagram of Fig. Q] is that of the solid line b. 
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FIG. 5: Scalar and pseudoscalar condensates vs. [id, for \x u = 200 MeV , T = (see Fig. QJ^. 
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FIG. 6: Phase diagram for chiral symmetry in the (/zb, T) plane for zero or small isospin chemical 
potential fj,i. For fj,i = (central line), the cross-over transition line starts from the point Pt = 
(0, 174) and ends at the point E = (220,85). The line between E and the point P M = (293,0) is 
the line for the first order transition with discontinuities in the (uu) and (dd) condensates. For 
fij = 30 MeV (side lines), the two cross-over transition lines start from the point Pt = (0,174) 
and end at the points E u = (190,85) and Ed = (250,85). The lines between E u and the point 
Pu = (263, 0) and between E& and the point P^ = (323, 0) are the lines for the first order transitions 
with discontinuities in the (uu) and (dd) condensates respectively. 
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FIG. 7: Phase diagram for chiral symmetry restoration in the plane {p^u-i^d) °f quark chemical 
potentials, at T = 60 MeV , which is below the temperature of the critical ending point (see Fig. 
EJj. Different regions are specified by the non vanishing of a given condensate, whereas the others 
are vanishing (p) or ~ m/A (\ u an d Xd)- Dashed lines are lines for the continuous vanishing of p 
or for cross-over phase transitions for \u or Xd? whereas solid lines are for discontinuous behaviors. 
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FIG. 8: Phase diagram for chiral symmetry restoration in the plane (fJ- u ,fJ^d) of quark chemical 
potentials, at T = 140 MeV, which is above the temperature of the critical ending point (see Fig. 
0). Different regions are specified by the non vanishing of a given condensate, whereas the others 
are vanishing (p) or ~ m/A (\ u an d Xd)- Dashed lines are lines for the continuous vanishing of p 
or for cross-over phase transitions for \u ° r Xd? whereas solid lines are for discontinuous behaviors. 
The solid line a refer to the path at ps = vs. pj > followed in Fig. \1(A 
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FIG. 9: Scalar and pseudoscalar condensates vs. fi^, for fi u = 200 MeV, T = 60 MeV (see Fig. 
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FIG. 10: Scalar and pseudoscalar condensates vs. fxj, for fiB = and T = 140 MeV . The figure 
is obtained following the path a in the phase diagram of Fig. 
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FIG. 11: Scalar and pseudoscalar condensates vs. fid, for \\i u = 200 MeV, T = 140 MeV (see Fig. 



0). 




FIG. 12: Critical value of the isospin chemical potential, beyond which a pseudoscalar condensate 
forms, vs. temperature, at zero baryon chemical potential. 



21 



— I — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — r 

fj, B = ; \±i = 200 MeV 

P 



0.5 



_i i i i I i i i i I i i i i Lli i i i I i i i i I i i i i 

50 100 150 200 250 300 

T {MeV) 



FIG. 13: Pion condensates vs. T , for hb = and jii = 200 MeV . The scalar condensates \u and 
Xd are of order ~ m/A. 



